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Abstract. From an analysis of the sky and velocity distributions of the high- velocity clouds 
(HVCs) we show that the majority of the HVCs has a common origin. We conclude that the 
HVCs surround the Galaxy, forming a metacloud of ~ 300 kpc in size and with a mass of ~ 
3x 10 9 Mq, and that they are the product of a powerful "superwind" (about 10 58 ergs), which 
, occurred in the Magellanic Clouds about 570 Myr ago as a consequence of the interaction 



of the Large and Small Magellanic Clouds. The HVCs might be magnetic bubbles of semi- 
ionized gas, blown from the Magellanic Clouds around 570 Myr ago, that circulate largely 



q . through the halo of the Galaxy as a stream or flow of gas. 

S_i ■ 

,!__> ■ On the basis of the connection found between the HVCs and the Magellanic Clouds, we 

C/2 . 

have constructed a theoretical model with the purpose of computing the orbits of a sample 
of test particles representing the HVCs, under the gravitational action of the Galaxy and 
the Magellanic Clouds. The orbits of the Large and Small Magellanic Clouds have been 
traced backwards in time to estimate the position and velocity of the Clouds at the time of 
the collision between the two Clouds, and to infer the initial conditions of the HVCs. The 
model can reproduce the main features of position and velocity distributions of the HVCs, 
like the overall structure and kinematics of the Magellanic Stream. The initial velocities of 
the HVCs were the result of velocities of expansion that permitted the escape of the HVCs 
from the Magellanic Clouds plus the systemic velocity of the Magellanic Clouds at the time 
of the collision. With these initial conditions, the Galactic gravitational potential induced 
differential rotations or shearing motions that elongated the cloud of HVCs in the orbital 
direction, forming the rear and front parts of the Magellanic stream. The population of 
HVCs is centered around the Magellanic Clouds. The eccentric position of the Sun within 
the cloud of HVCs explains the asymmetries between the sky distributions of the HVCs of 
the northern Galactic hemisphere and those of the southern Galactic hemisphere. 

In the light of the model we analyze the effects that the passage of the HVC flow through 
the Galactic disk has produced on the interstellar medium. The effects of the HVC flow can 
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account for many observational details such as the Galactic warp, HI shells and supershells 
in the gaseous layer of the outer parts of the Milky Way. The Galactic disk was target of 
numerous impacts of HVCs in the course of the last 400 Myr, accumulating mass at the 
average rate of approximately 0.6 Mq per year. The events of this period may be regarded 
as landmarks in the evolutionary history of the Milky Way. 

Key words. ISM: clouds - Magellanic Clouds - Galaxy: structure - Galaxy: halo - Galaxy: 
evolution - galaxies: interactions 

1. Introduction 

The existence of HI high-velocity clouds (HVCs) at high Galactic latitudes has intrigued as- 
tronomers from the time of their discovery in 1963. The early observations induced the belief 
that all HVCs had negative radial velocities, and that these clouds could be intergalactic gas 
falling onto the Galaxy. The recognition of violent activity in the Galactic disk and the discovery 
of large complexes of high positive velocity clouds and of the Magellanic Stream opened a whole 
gamut of possibilities for the nature and origin of the HVCs such as "Galactic fountain" clouds, 
material stripped from the Magellanic Clouds and extragalactic objects. The lack of adequate 
distance indicators makes it difficult to discriminate among the different alternatives. 

The first interpretations of the HVCs were discussed by Oort J 19661 119691 11970> . Wakker 
& van Woerden (I1997> have given a recent review of the HVC phenomenon. Poppel J1997I 
reviewed the possible role played by the HVCs in the local interstellar medium. Although a large 
amount of observational data has now been accumulated, such as very sensitive HI 21 cm surveys 
(Morras et al. BOOOl de Heij et al. l2002all2002bl Lockman et al. l2002l Putman et al. l2002l Putman 
et al. 2003ai, the detection of molecules (Richter et al. 20011 and Ho- emission in HVCs (Tufte 
et al l 19981 Putman et al. l2003b> and the distance constraints to some HVCs (Wakker l2"0TTlV as 
well as a large number of theoretical studies that have been dedicated to this topic (eg. Quilis & 
Moore l200Tl Espresate et al. l2002l Sternberg et al. l2002l Maloney & Putman 120031 . the HVCs 
remain enigmatic. 

The theories for the origin and nature of the high velocity clouds can be divided into three 
categories: intergalactic (eg. Bajaja et al. 119871 Blitz et al. 119991 Braun & Burton |2000> . cir- 
cumgalactic (eg. Kerr & Sullivan III 1969 Hulsbosch & Oort 11973l and Galactic (eg. Bregman 
1980, Verschuur 1993). On the basis of an analysis of the sky and velocity distributions of the 
HVCs we will demonstrate in this paper that the majority of these clouds constitutes a circum- 
galactic flux or stream of clouds related to the Magellanic Stream. The idea that the HVCs may 
be fragments of Magellanic material precipitating toward the galactic disk was first suggested by 
Giovanelli ( 1981 1 and Mirabel ( 1981 1. The main theories on the origin of the Magellanic Stream 
itself consider that the stream consists of gas extracted from the Magellanic Clouds. The pro- 
Send offprint requests to: C. A. Olano 
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posed mechanisms that have been studied in detail are tidal forces of the Galaxy, friction forces 
of the gaseous galactic halo and a collision between the Large and Small Magellanic Clouds 
(Gardiner et al. lT9951 Heller & Rohlfs[T333| Moore & Davis Gardiner & No guchi [19961 . 
The tidal forces affect equally both gas and stars of the Magellanic Clouds. However there are 
no stars in the Stream, indicating that tidal interaction could not be the primary cause of its for- 
mation. On the other hand, the weakness of the stripping drag hypothesis is that the gas density 
of the Galactic halo should be very low. Hence, to account for the formation of the Magellanic 
Stream and the rest of the HVCs, we will here take a different point of view that includes only 
the collision between the Large and Small Magellanic Clouds as the mechanism that triggered 
the process. 

2. A simple interpretation of the sky distribution and kinematics of the high 
velocity clouds 

In general the observations of HVCs provide the HVC positions projected on the sky and their 
radial velocities. Hence, we do not know the distances and tangential velocities of the HVCs. 
Even with this limitation, if the HVCs are a homogeneous population with a common origin, in 
principle it is possible to find the space distribution and the three-dimensional velocity distribu- 
tion of the HVC population from the analysis of the observed distributions of the sky positions 
and radial velocities of the HVCs. This is a typical case of an inverse problem in astronomy 
(Craig & Brown ll986l Merritt lT996l Saha ll998l . 

Nowadays a vast amount of observational material is available on HVCs in the form of sur- 
veys whose combination covers the whole sky, and that therefore provide complete and reliable 
statistical information. Wakker ( 2002 1 compiled a catalog of HVCs with 1 1000 entries, using the 
surveys of Hulsbosch & Wakker (JT988), Bajaja et al. (IT9851 for -23° < 6 < -18° and Morras et 
al. J2000> for 6 < -23°. Each entry contains the celestial position coordinates, I and b, the radial 
velocity p and other parameters that characterize a detected high-velocity profile component (i.e. 
a piece of cloud intercepted by the radio telescope beam). An HVC can be defined by a set of 
points {(, b, p), or cloud elements, that satisfy criteria of continuity in the (€, b, p)-space. For the 
study of HVC distributions, we will represent all the HVC components as separate points (f, b, 
p); although each point can be associated with a subset of points (defining an HVC or complex of 
HVCs) that are not necessarily independent. In statistical considerations one should take this into 
account (see Sect. 4). We will represent the distributions of HVCs in the conventional form as 
orthogonal projections of points in the ({, b, p)-space upon the planes l-b, l-p and b-p, although 
a three-dimensional representation is more elegant and compact. 

On the basis of the data of Wakker's catalog, we display the distributions of the positions and 
radial velocities of the HVCs in Figs. 1 1 -al II -bl and ll-cl For purposes of comparison, in 
Figs. l2-al 12 -bl and 12-cl we represent the HVC distributions corresponding to a second set of data 
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Fig. 1-a. Aitoff projection of the entire sky distribution of the HVCs in Galactic coordinates, 
based on measurements of the 21-cm emission line of HI (adapted from Bajaja et al. 119851 
Hulsbosch & Wakker 1988 and Morras et al. 2000 1. Plus (or darker regions) and minus signs 
mark respectively the positions of the HI lines with positive and negative high velocities in the 
LSR frame. 

with an improvement in spatial resolution taken from Putman et al. ( 2002 1 and Giovanelli ( 1980) 
for 5 > -15°. 

The enclosing curve of the longitude-velocity distribution (Figs. II -bl and 12-bl is approx- 
imately sinusoidal and that of the latitude-velocity distribution (Figs. II -el and 12-cl follows a 
cosine law. In consequence, the radial velocity distribution of the HVCs can be represented by 
the following formula: 

p = -V sin(€)cos(b), (1) 

where V - 450kms~'. This distribution law admits a simple interpretation (or "inversion"), 
namely: it reflects a roughly uniform flow of HVCs that comes approximately from galactic 
co-ordinates £q = 90° and bo = 0° and that encloses the whole Galaxy. If Vq and Vf are the 
velocities of the Sun and of the HVC flow with respect to the Galactic center respectively, the 
radial velocities of the HVCs with respect to the LSR are given by 

p = -Vq sin({) cos(b) - V f cos(t - £ Q ) cos(b - b ). (2) 

Substituting the values of (q and bo into Eq. (0, we obtain Eq. Q with V = Vq - Vf. The velocity 
of the HVC flow is of the order of the Solar velocity (220 km s~ 1 ), and has the same direction, but 
its sense is reversed. In other words, the Sun is immersed in the flow of HVCs, moving counter 
stream. This flow of clouds seems to circulate largely in the halo or in the outskirts of the Galaxy. 
In the next section, we will find a clue that makes it possible to connect the flow of HVCs with 
the Magellanic Clouds and the Magellanic Stream. 
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Fig. 1-b. Velocity of high-velocity HI lines with respect to the LSR plotted against Galactic lon- 
gitude, from the collection of data used in Fig. Il-al The sinusoidal curve represents the Vlsr(^) 
pattern for b — 0° corresponding to a hypothetical flow of HVCs, circulating in directions parallel 
to { = 90° - 270°. 

3. Model for the origin and dynamical evolution of the high velocity clouds 

In this section we will try to determinate the probable origin of the HVC flow described in Sect. 
2 and the initial conditions of the phenomenon; with this we will construct a model to study the 
trajectories of the HVCs throughout the halo and the disk of the Galaxy. Our first objective is to 
show that the Magellanic Clouds have played a central role in the origin of the HVCs. For this 
purpose, we need to calculate the orbital plane of the Magellanic Clouds. We adopt throughout 
this paper a Cartesian coordinate system (X, Y, Z) with the origin at the Galactic center, the X-axis 
pointing in the direction of the Sun's Galactic rotation, the Y-axis pointing in the direction from 
the Galactic center to the Sun, and the Z-axis pointing toward the Galactic north pole. Since the 
Magellanic Clouds are subject to a central force due to the Galactic halo (we will ignore the forces 
of the galactic disk), their orbital plane can be characterized by the vector normal to this plane, 
namely N = r xr, where r - (X, Y, Z) is the position and r — (X, Y, Z) is the velocity vector of the 
center of mass of the combined Clouds at the present epoch in the Galactocentric rest frame. We 
use the observational parameters of the LMC, which are better determined, as representative of 
the Magellanic Clouds as a whole. We adopt the following parameters for the LMC: i = 280°. 5 
and b = -32°. 9 for the Galactic coordinates, 50. 1 kpc for the distance from the Sun, 262.2 km s 
for the heliocentric systematic velocity (van der Marel et al. 20021, fi a = 0.00486" yr _I and 
fig = 0.00034" yr~' for the mean proper motion (see Table 1 of van der Marel et al. 12002b . 
To correct for the reflex motion of the Sun and to obtain the positions and velocities in the 
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Fig. 1-c. Velocity of high-velocity HI lines with respect to the LSR plotted against Galactic lat- 
itude, from the collection of data used in Fig. Il-al The two half-cosine curves represent the 
Vlsr(£) pattern for I = 90° and 270° corresponding to a hypothetical flow of HVCs, circulating 
in directions parallel to I = 90° - 270°. 



Galactocentric frame of rest, we use the IAU values Rq = 8.5 kpc and Vq = 220 km s _1 and the 
"basic" solar motion of 16.5 km s _I towards the direction of i — 53° and b = 25° (Allen IT9731 . 
With these adopted values, we calculate the position and velocity vectors of the LMC 



r lmc = (-41.30, 1.29, -27.43) kpc (3) 
and 



W = (-218,66, 187) km s" 1 . (4) 

The resulting plane of the orbit of the LMC is represented in Fig. [3] The angle between the 
orbital plane and the plane perpendicular to the Galactic disk is only of the order of 10°. i.e the 
plane of the LMC orbit is almost perpendicular to the Galactic plane. The Galactic longitude of 
the nodal line is nearly 82°, and thus the orbital plane is almost coincident with the X-Z plane (cf. 
Fig. 1 of Murai & Fujimoto 1980 1. Therefore, the three-dimensional velocity vectors representing 
the HVC flow described in Sect. 2 are nearly parallel to the plane of the orbit of the Magellanic 
Clouds. This indicates a dynamical relationship between the HVCs and the Magellanic Clouds. 
Our working hypothesis will be that the HVCs were ejected from the Magellanic Clouds as a 
result of internal processes in the Magellanic Clouds, activated likely by the collision of the 
LMC and SMC. 
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Fig. 2-a. Same as Fig. Il-al but with data adapted from Giovanelli ( 1981 1 and Putman (12002> . 
The region centered at I — 120° and b = 30° is not covered by the observations of the two 
surveys used here, and is left empty. 
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Fig. 2-b. Same as Fig. Il-bl but from the collection of data used in Fig. 12-al 

3.1. Basic equations of the model 

The basic problem we should consider here is the motion of a test particle representing an HVC 
in the gravitational fields of the Galaxy and the Magellanic Clouds, and friction forces due to the 
gaseous disk of the Galaxy. 




In our computations, we use the spherical gravitational potential of the halo iff = - V 2 In r for 
the Galaxy (Murai & Fujimoto 1980, Lin & Lynden-Bell 1982), which gives a flat rotation curve 
with a constant circular velocity, V c = 220kms~ 1 , so that the gravitational force of the Galaxy 
exerted on a particle of unit mass is -V 2 4. For the Magellanic Clouds, we use the gravitational 
potential of the LMC and ignore that of the SMC. We consider the SMC as another test particle 
(see Sect.3.2). We assume that the LMC has a Plummer-type potential with an effective radius 
K = 3kpc (Murai & Fujimoto 1 19801 . giving a gravitational force per unit mass of - ™^^ j2 , 
where r'{- r — j"lmc) denotes the position of the test particle with respect to the LMC, and Mlmc 
is the total mass of the LMC estimated at 0.87 x 10 10 M Q (van der Marel et al. 2002 1. 

The friction force per unit mass that acts on an HVC can be expressed as F A = -£^-V r V r , 
where v r is the relative velocity between the HVC and the gaseous medium in which the cloud 
moves, n is the density of the gaseous medium at the position of the HVC, S is the head cross- 
section of the HVC, m is the HVC mass, and s is a dimensionless quantity 1. A measure 
of the ratio j is the column density Nh of the HVC being studied, an observable quantity. We 
adopt a mean Nh of 7.2 x 10 18 atoms crrT 2 for the HVCs. In our simplified model, the Galactic 
halo is empty of gas, thus n depends only upon the density distribution of the gaseous disk of 
the Galaxy. The density n may be conveniently expressed in cylindrical coordinates as n{R, z) = 
exp[- J^], where cr(R) and h(R) are the surface density of interstellar HI and the scale 
height of the thickness of the HI gas layer at R (= ^x 2 + y 2 ), respectively. We denote by x, y and 
Z the Cartesian components of the position vector r of the test particle or HVC. From a fit of 
Wouterloot et al.'s (fT990l Table 1 we derive cr(R) = 1.0875 x 10 21 exp[-( ^g^ ) 2 ] atoms cnr 2 
and h(R) = 105.8 exp[( 4 "^ ? ) 2 ] pc with R in kpc. Assuming that the gas of the Galactic disk 
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Fig. 3. Locations of the LMC, the SMC and the Sun in Galactocentric coordinates (X, Y,Z). The 
past orbit of the LMC is represented by the solid curve. The position angles of the orbital plane 
are indicated. 



rotates circularly with constant velocity V c , the velocity vector of the gaseous medium at the 
position of the HVC is (-^y, — ^jc,0) and hence v r = (x- ^jy^ + ~jf x >£)- Under these conditions, 
the equations of motion of a test particle are 

.. y]r GM LMC (r - r LMC ) _ 

r+ ^ + (\r-r LMC t + K^ - F « (r > r) = °- (5) 

To solve this system of equations, we need first the solution of the equations of motion of the 
LMC in the gravitational field of the Galaxy, namely: 

V^r LMC 

f-LMC + — 2 = 0. (6) 

r LMC 
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3.2. The collision between the Large and Small Magellanic Clouds 

Several authors agree that the Magellanic clouds had a close encounter (separation less than 3 
kpc) between 200 and 500 Myr ago (Murai & Fuiimoto lT9"gUl Gardiner et al. [13331 Heller & 
Rohlfs lT994l Moore & Davis ll994l Gardiner & Noguchi [T996t . The morphological peculiarities 
of the Clouds such as the widely scattered distribution of HII regions around the optical bar of 
the LMC and the severe fragmentation of the SMC could be understood in terms of this collision. 

The encounter time f e of the Clouds is an important parameter in our model, and we want to 
calculate its value congruently with the other parameters adopted above. The orbit of the LMC 
can be computed straightforwardly by means of the numerical integration of Eq. l|6} with the ini- 
tial conditions given by Eqs. (0 and (0} (see Fig. [3}- For simplicity we assume that the orbit of the 
LMC is not essentially altered by the presence of the SMC because of its lesser mass. Knowing 
the present position and velocity of the center of mass of the SMC, we compute the orbit of the 
SMC by using Eq. (0 with Fa — and the solution of the orbit of the LMC. Thus we obtain the 
time of the closest approach of both orbits. However, the position and velocity of the SMC are 
not known with great precision. Therefore, in the equations we only introduce the radial compo- 
nent of the SMC velocity, based on the heliocentric radial velocity of 148 kms~' for the SMC 
that was accurately measured with optical and radio-astronomical methods (Hardy et al. U989> . 
and leave as unknowns the proper motions and the distance, which were used for determining the 
tangential components of the velocity. Solving Eq. (0 with the help of Eq. we obtain tsmc 
as a function of fj.„, fig, the distance d and time t. By requiring that rsMc(^a,Ms, d, t) = riMc(t)> 
we have three equations with three unknowns for each chosen past time; their solutions give us 
fj. a (t), /us(t) and d(t). Choosing the set of theoretical parameters compatible with the observational 
ones, we found that fi a , fig and d should be close to 0.00432" yr _1 , -0.00108" yr~' (cf. Kroupa & 
Bastian 1997) and 60.3 kpc (cf. Westerlund 1990 1 respectively, and that the collision between the 
Clouds occurred 570 +350 Myr ago. The uncertainty of +350 Myr in the value of t e was derived 
from that of the present velocity of the LMC (Eq. 4), Aruvic ~ (±23, ±36, ±35)kms~ 1 (see Eq. 
55 of van der Marel et al. 2002 ), with the assumption that the rest of the adopted parameters like 
the distances to the Large and Small Magellanic Clouds and the gravitational potential used for 
the halo of the Galaxy are almost exact. The orbits of the Clouds are displayed in Figs. 14-al and 

The velocity of the LMC at the time of the collision f e = -570 Myr was 

/•LMc(fe) = (-100, -12, -141) km s-\ (7) 

while that of the SMC was rsMcfe) = (-100, -167, -155)kms _1 . Thus, during the encounter, 
the velocity of the SMC relative to the LMC was of order 156kms~'; i.e. very close to the 
escape velocity from the LMC (v e = -^ 2G ^ LMC - 160 km s -1 ). To explain this, we propose that 
the Magellanic Clouds formed a binary system in the past and that it was broken by a sudden 
mass loss of the system taking place in the epoch of the collision. The collision compressed 
the interstellar gas of the Clouds, inducing a burst of star formation. The massive stars of this 
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Fig. 4-a. Projection of the past orbits of the LMC and the SMC on the Galactocentric X-Y plane. 
The LMC and SMC positions are marked by crosses each 100 Myr. 



Z(kpc 




Fig. 4-b. Projection of the past orbits of the LMC and the SMC on the Galactocentric X-Z plane. 
The LMC and SMC positions are marked by crosses each 100 Myr. 



new stellar generation expelled the surrounding interstellar matter in a few millions of years 
at velocities greater than the escape velocity of the system, through supernova explosions and 
strong stellar winds. A part of the material ejected from the Clouds at that time formed the 
HVCs. The HVCs left likely the Magellanic Clouds as magnetized clouds of plasma (eg. Olano 
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119991 ). Magnetic fields in HVCs have been detected by measuring Zeeman splitting of the 21- 
cm HI line (Kazes et. al. 1991 1. This magnetic self-confinement of the HVCs could explain 
the HVC longevity as discrete entities. This idea is coherent with the fact that the HVCs are 
highly ionized (Sembach et al. [T9^9l Lehner et al. l200Tl Sembach et al. 120051 . The HVCs are 
still generally thought of as neutral entities. However, at an early stage they were probably fully 
ionized. Subsequently, the evolution of their physical conditions permitted the recombination 
of part of the protons and electrons. On the other hand, the explanations generally given for 
the cloud confinement are that they are either bound by gravitationally dominant dark matter or 
confined by an external pressure provided by a hot gaseous Galactic halo (Espresate et al. 120021 
Stanimirovic et al. 120021 Sternberg et al. 2002 1. 



3.3. Initial conditions of the high velocity clouds 

We conjecture that the HVCs were ejected from the Magellanic Clouds as the consequence of 
very energetic interactions of a large number of stars born during the collision of the Clouds 
and the cloudy interstellar medium of the Clouds. The relevant initial conditions, namely the 
positions, velocities and times of the ejections of the HVCs, should therefore be characterized 
by distribution functions. For simplicity we assume that the HVCs were ejected from a relatively 
small region of the Clouds, at the central point of the collision, during a short period of time, i.e. 
at f e =-570 Myr. The position vector of the central point of the collision (see Sect. 3.2 and Figs. 
14-aland 14-bb is then 



/-LMc(0 = (86.39, -15.53, -15.68) kpc, (8) 

and the initial position vector for all HVCs is ro = fLMc(fe). For the initial systematic or 
mean velocity of the HVCs we adopt the velocity of the mass center of the LMC at the 
time f e , which is given by Eq. @. The errors in Eqs. @ and (8), i.e. in the position of the 
LMC in phase space at t e , are Ar LM c(fe) ~ (±3.87, +0.72, +1.03)x | Ar LM c I kms~' and 
A/"LMc(fe) ~ (±0.79, ±0.20, ±1.21)x | A/lmc I kpc, obtained by propagation of the errors in 
/•lmc an d t e (Sect. 3.2). 

Thus the initial total velocity of an HVC can be written as 

ro = r LMC (t e ) + v p , (9) 

where v p is the initial peculiar velocity of the HVC. The number of HVCs in the initial velocity 
volume dv p at v p is dN = f(v p )H(v p - v e )dv p at time f e . We assume that the distribution function 
/(v p ) follows Schwarzschild's ellipsoidal law: 

W = n ^3/2v v v exp -i((^) 2 + (^) 2 + (^) 2 )), (10) 

{Z7Ty^LiL2^3 2 2-1 2,2 2-3 

where Si, £2 and S3 are the three velocity dispersions in the directions of the three co-ordinate 
axes and i p , y p and z p are the components of the peculiar velocity vector v p . The Heavyside 
function H, which is 1 for v p - v e > and otherwise, takes into account that only the clouds 
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with velocities greater than the escape velocity v e became HVCs. In general form, the func- 
tion of initial distribution of the HVCs may be written as F(ro, ro, t) = <5(ro - fLMc(fe))<5(f - 
f e)/('"o _ r LM c(t e ))H(ro - r LM c(f e ) - v e ), where 6 is the standard delta function. The total num- 
ber of HVCs (or the number of test particles used in the simulation) can be expressed by 
Nhvc = C° F(ro, ro, f)dr$drodt, which makes it possible to determine N of Eq. (II Oi . 

4. Results of the model and comparison with the observations 

To reproduce the main features of the position and velocity distributions of the HVCs we cal- 
culated the orbits of about 850 test particles, which represent the HVCs. Distributing these test 
particles in phase space according to the distribution functions defined in Sect. 3.3, we can as- 
sign an initial position and velocity to each test particle, and calculate its present position and 
velocity by means of the equations of motion given in Sect. 3.1. The three dispersions of the 
velocity ellipsoid in Eq. (II Oi are the free parameters of the model. Their values can be estimated 
by comparing the predictions of the model with the observations. We found that the case par- 
ticular of a three-dimensional Gaussian distribution is compatible with the observations, namely 
Ei = 1.2 — E3 = E = 120 km s _1 . Thus, E is the only free parameter of the model. 

To consider the possibility of non-isotropic expulsion of material, we tried other initial el- 
lipsoidal distributions with different values and spatial orientations of their velocity dispersions 
with respect to the co-ordinate axes, but none of them could produce better results than those of 
a Gaussian distribution. This simple distribution reproduces fairly well the kinematics and sky 
distributions of the Magellanic stream, the "classic" stream and the leading stream (Putman et al. 
2003a), as well as those of the rest of the HVC population. Certainly, quantitative comparisons 
between theoretical and observed distributions in the way developed by Saha ( 1998 1 would be 
desirable. However, it is beyond the scope of the present work. A proof of the robustness of our 
method is that the results of the simulations are rather insensitive to the uncertainties of the time 
of the collision, Af e , and the velocity, AruvicfeX an d position, ArLMc(fe), of the centroid of the 
initial Gaussian distribution. 

The results of the model are shown in Figs. 15 -al 15 -b I 15-cl 15 -dl and 15 -el Comparing Fig. 
15-al with Figs. II -al and 12 -al Fig. 15 -bl with Figs. II -bl and 12 -b I and Fig. 15-cl with Figs. Il-cl 
and 12-cl we see that the theoretical distributions are in good agreement with the corresponding 
observed ones. An exception is a group of HVCs that lies in a region centered at I — 180°, 
b - -20°, whose velocities are not explained by the model. However, we should remember the 
simplifications and uncertainties of the observational parameters introduced into the model, as 
well as the probable existence of other sources of HVCs, apart from the Magellanic Clouds. 
Studies of HVCs in external galaxies (see Schulman et al. 1996, Wakker & Woerden 119971 
Jimenez- Vicente & Battaner 2000, Miller et. al. 2001) will enable us to make a comparison 
with the HVC system associated with our own Galaxy and to decide the relevance of the other 
possible sources of HVCs, such as Galactic fountains (Bregman 1980 Houck & Bregman 1990). 
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The part of the HVC population of Magellanic origin having LSR radial velocities between - 
100 and +100 kms -1 and lying at low Galactic latitudes (see Figs. 15-bl and 15-cl should be 
indistinguishable from the interstellar clouds of the Galactic disk that are parts of supernova 
shells or that participate simply in the Galactic rotation. 

There is a clear difference between the two Galactic hemispheres. In the northern hemisphere 
the population of HVCs is dominated by a few extended, complex structures (the A, C and M 
objects). When comparing the observed HVC distributions with the theoretical ones, we should 
remember this. An HVC sampled at regular angular intervals is characterized by a set of neigh- 
boring points in the (€, b, p) space whose number increases proportionally with the solid angle 
subtended by the cloud in the sky (see Sect. 2). The few northern complexes of HVCs cover a 
large area of the sky. Hence their densities of points in the observed distributions relative to the 
densities of points in the corresponding theoretical distributions should be greater than the rela- 
tive densities of points corresponding to the southern HVCs. The results of our simulation agree 
with this picture (cf. Figs. II -al and 15 -al and Figs II -el and 15-cl . According to our model, the 
explanation of these asymmetries between the northern and southern hemispheres is that in the 
northern Galactic hemisphere there are far fewer HVCs and they lie at much smaller distances 
(see Fig. 15-el . Therefore, they subtend far larger angles on the sky. The importance of the model 
is that it provides a probable velocity field and space distribution of the HVCs in the three di- 
mensions (Figs. l5~^dl and 15-el . This information cannot be obtained from the observations so far. 

The group of HVCs spreads throughout a volume of approximately cubic 300 kpc. The mean 
velocity of the HVC flow is similar to that of the Magellanic Clouds, since the HVCs started their 
trajectories with the mean or systematic velocity of the Clouds. To illustrate this effect further, 
we apply the model to a set of test particles with initial peculiar velocities of the same magnitude 
| v p |=cont, but with different directions. At the present time, these test particles should form 
a stream that moves approximately in the same direction and sense as the Magellanic Clouds 
(Fig. 16-al . In addition, this cloud of particles is elongated in the direction of the velocity of 
rotation around the Galactic center, an effect due to the differential rotation and well known 
in the evolution of expanding Galactic shells (e.g. Olano 1982; Palous, Franco & Tenorio-Tagle 
1990] Olano|2001 1. Notice that the degree of concentration of the particles around the Magellanic 
orbit increases as the initial peculiar velocity decreases (cf. Figs. 16 -al and 16-bl l. Also, notice that 
the lower the peculiar velocities of the particles, the greater the number density of particles (or 
HVCs), according to the initial distribution function Eq. JlOl i. Both these effects contribute to 
create in the sky the appearance of the Magellanic Stream (see Fig. 15-al . 

From the distances to the HVCs predicted by the model and the density column data of 
Wakker's ( 2002 ) list we estimated the total mass of the HVCs by means of the formula M x (Mq) = 
6.085 x lO-^Q^i Af H .), where N H , is the sum of the column densities in atoms cm of 
all high-velocity profile components detected with a telescope beam of ~ 0°.5, and cfi is the 
theoretical mean squared distance of the HVCs from the Sun, expressed in kpc 2 . The results are 
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Fig. 5-a. Simulated celestial distribution of the HVCs: projection of the positions of the test parti- 
cles at present moment on the sky plane. The plus and minus symbols have the same meanings as 
in Fig. Il-al The projections of the past orbit of the LMC (solid curve) and the present positions 
of the LMC and SMC (large and small projected circles) are shown. The cross at the end of the 
projected orbit curve is the projected position of the encounter between the LMC and the SMC. 
Compare this theoretical distribution with the observed ones (Figs. II -al and 12-at . 
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Fig. 5-b. Simulated velocity-longitude relation for the HVCs: radial LSR velocities of the test 
particles at the present moment as a function of their respective Galactic longitudes. The solid 
curve has the same meaning as in Fig. II -bl Compare this theoretical relation with the observed 
ones (Figs. II -bl and 12-bi . 
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Fig. 5-c. Simulated velocity-latitude relation for the HVCs: present radial LSR velocities of the 
test particles as a function of their respective Galactic latitudes. The solid curve has the same 
meaning as in Fig. Il-cl Compare this theoretical relation with the observed ones (Figs. II -el and 
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X? =1 N Hi = 1. 87 xlO 23 atoms enr 2 and yjdf = 105 kpc, hence the total HI mass of the HVCs is ~ 
1 .2 x 10 9 Mq. Assuming that the HVCs have as much ionized as neutral hydrogen (Sembach et al. 
2002i, and adopting a factor 1.3 to include H e , we obtain M t = 3.1x10 9 Mq. Thus the Magellanic 
Clouds lost about 25 per cent of their original mass in the form of HVCs. This is consistent 
with the idea that the Clouds constituted a binary system before the last collision, whereas the 
important mass loss induced by this collision transformed the Clouds into an unbound system. 

To study the history of the Clouds as a binary system we should calculate the relative motion 
of the SMC and LMC before the collision. This can be realized by means of the equation of mo- 
tion resulting from the difference of Eq. (0 and Eq. (|6j, the boundary conditions at the moment 
of the collision, and a LMC mass greater than the present one (see Sect. 3.1). However, for this 
purpose we merely quote the pioneering study on this topic (Murai & Fujimoto 1980 1. According 
to it the LMC and the SMC formed a stable binary system during 10 10 yr (see also Yoshizawa & 
Noguchi l2003l . 

Another interesting result of the model is the estimate of the total energy associated with 
the production of the HVCs. The total initial kinetic energy of all HVCs (jM t Vp) is about 1.8 x 
10 57 ergs, which implies processes that generated at least a total of 10 58 ergs. This huge amount 
of energy should have been supplied by supernovae and winds from massive stars in a starburst 
originated by the collision of the Clouds. The collective action of supernovae and stellar winds 
of a typical starburst can drive a "superwind" that may eject 10 7 - 1O 9 M0 from a galaxy with 
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Fig. 5-d. Simulated spatial distribution and velocity field in the Galactocentric rest frame for 
the HVCs: Present positions and velocities of the test particles. The test particles with initial 
velocities | v p |< 220 km s _I have been omitted in the plot for clarity. The present positions of the 
LMC and the SMC (circles), the past orbit of the LMC (solid curve) and the point of encounter 
between the LMC and the SMC (circle) are represented. 



a mechanical energy of 10 57 - 10 59 ergs during its lifetime estimated at 10 7 yr (Heckman et al. 
1989 1. There are various stellar generations in the Magellanic Clouds (Westerlund 1990 1. There 
is nothing against the idea that the formation of the younger generations, with ages lower than 
600 Myr, was initiated by the powerful starburst associated with the collision of the Clouds about 
570 Myr ago. 
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Fig. 5-e. Same as Fig. 15 -dl but projected on the Galactocentric X-Z plane. 
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5. Interaction of the high velocity clouds with the Galactic disk and 
generation of the Galactic warp 

Most of the HVCs have not been affected by the friction force of the interstellar medium, since 
their orbits have not intersected with the gas layer of the Galactic disk. Here, we analyze the 
important group of HVCs that was strongly decelerated by the friction force of the Galactic layer 
and that had a strong influence on the Galactic disk. 

The HVCs that penetrated the dense regions of the gaseous disk were trapped in the Galactic 
plane, and are participating in the Galactic rotation now (see Figs. 17 -al and 17-bl i. In Fig. [8] 
we show the positions of entry of the test particles in the gaseous layer of the Galactic disk and 
their present positions, resulting from the braking process. According to our model, during the 
last 400 Myr the Galactic disk has been profusely hit by HVCs (see Fig. [8]and Fig. 19-at . 

The detailed description of the interaction of HVCs with the Galactic disk is beyond the scope 
of the present study. Tenorio-Tagle ( 1980*l ll98Tt analyzed the physics of the collisions of HVCs 
with the Galactic disk. It is clear that the impact of HVCs upon the Galactic disk should be an 
important source for producing shells and supershells. This mechanism of creation of large-scale 
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Fig. 6-a. Spatial distribution and velocity field of the test particles with an initial isotropic veloc- 
ity distribution characterized by | v p |= 170 km s _1 . 

structures in the interstellar medium can account for many of the observational details of the large 
HI shells lying beyond the solar circle (Mirabel 1982 Heiles 1984 1. In connection with this, it is 
interesting to compare the theoretical distribution of the impacts of test particles on the Galactic 
disk (Fig. [HJ, with the location of observed shells and supershells of HI (McClure-Griffiths et al. 
2002 1. Evidence for interactions of HVCs with the gas in the Galactic plane was found by several 
authors (see Tripp et al. 2003 and the reviews of Tenorio-Tagle & Bodenheimer 1988 and Poppel 

nw7i . 

In the next paragraphs of this section we will show that the flow of HVCs onto the Galactic 
disk over a period of 400 Myr has caused a major disturbance in the Galaxy. As a consequence of 
the fall of clouds, a burst of star formation involving the whole Galactic disk should have started 
400 Myr ago and is still going on (Barry 1988 Noh & Scalo 1990 1. There are observational 
details of the vertical structure of the Galactic disk, such as peculiar departures from the midplane 
or huge depressions delineated by young star-gas supercomplexes (Alfaro et al. 1991 1, that may 
be explained in terms of HVC-disk collisions that induced star formation (Cabrera-Cano et al. 
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Fig. 6-b. Spatial distribution and velocity field of the test particles with an initial isotropic veloc- 
ity distribution characterized by | v p |= 270kms~'. 

1995 1. Further predictions provided by our model are the times, velocities and angles of entry 
of the HVCs in the Galactic layer (see Figs. 19 -al and 19-bt . The angle of entry (or incidence) 
of an HVC is defined as the angle of the direction of motion of an HVC with respect to the 
vertical (Z-axis). The angle of entry is an important parameter for the hydrodynamic simulations 
of the collision process (Comeron & Torra ll99"Tl Franco et al. 11988b . The times, velocities and 
angles of entry correspond to the point at which the friction coefficient reaches the threshold 
value s— = lO^pc' 1 (see Sect. 3.1). Figs 19 -al and 19-bl show interesting differences between 
the times and angles of entry on the right side of the Galactic disk (positive X-axis) and those on 
the left side. The HVCs arrived first at the right side of the Galactic disk, while there was a delay 
of ~ 150 Myr for those arriving at the left side of the disk. On the right side of the disk the HVCs 
impinged on the disk from below with angles of entry ~ 80°, i.e. almost a grazing incidence. In 
contrast, on the left side of disk the angles of entry are lower; and although most of the HVCs 
entered from below too, there were some HVCs that entered from above. 
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Fig. 7-a. Projection, on the Galactocentric X-Y plane, of the orbits of some of the test particles 
that interact with the Galactic disk. The filled circles indicate the positions of entry into the 
Galactic layer. The crosses mark the positions of the particle at intervals of 30 Myr. The times 
of entry of this subset of test particles lie between -220 and -120 Myr. The position of the Sun 
(X Q , Y q ) is (0, 8.5)kpc. 

In the following we examine whether there is a causal relationship between the HVC-disk 
interaction and the generation of the Galactic warp. Galactic warps are very complex phenom- 
ena (see Lopez-Corredoira et al. 2002 Garcfa-Ruiz et al. 2002 and references therein). We will 
analyze only a few aspects of the problem in the light of our model. During the last 450 Myr, 
the flow of HVCs has been exerting a pressure on the gaseous disk which could have caused the 
distortion of the disk of the Galaxy. Since about 7 per cent of the test particles of the simulation 
impacts within the galactocentric radius R = 25 kpc, the mass in the form of HVCs accumulated 
by the disk is 0.07 x M t ~ 2.2 10 8 M©. Therefore, the rate of mass accretion per surface unit is 
y = 600M©kpc~ 2 Myr~', taking into account that the time during which the flow acted on the 
disk at a fixed point in the (X, Y,Z) frame was nearly 200 Myr (see Fig. l9-at . Let us consider a 
cylinder of unit cross-section with generators parallel to the Z-axis and height h{R) as a volume 
element of the gaseous disk of the Galaxy at R. The accretion process acting on the volume ele- 
ment generates a force in Z per unit of mass that can be written as F z = ^^vcos(/), where v, i 
are the mean velocity and angle of entry, and cr(R) is the surface density of the gaseous disk at R 
(see Sect. 3.1). 

To describe the gas motion of the volume element under the effects of the additional force 
associated with the mass accretion, we use cylindrical coordinates, {R, 9, Z), representing the 
galactocentric radius, azimuth and vertical distance from the plane b = 0° of the mass center of 
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Fig. 7-b. Projection of the orbits of the same test particles as those of Fig. 17 -al but on the 
Galactocentric X-Z plane. 
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Fig. 8. Same as Fig. 17 -al but for the total set of test particles that interact with the Galactic 
layer. The filled and open circles indicate respectively the positions of entry into the Galactic gas 
layer, and the present positions of the impacted regions, resulting from the interaction of the test 
particles with a rotating disk. We have omitted the projection of their complete orbits for clarity. 
The numbers inside the circles refer to the ages of impacted regions in units of 10 Myr, defining 
age as the time elapsed from the entry of the particle into the Galactic layer up to the present. 
The solar circle is indicated. 
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Fig. 9-a. Correlation between the times of entry of the test particles into the Galactic layer and 
the positions of the particles projected on the X-axis. The solid lines represent the position of the 
front edge (lower line) and the back edge of the accretion front as a function of time. Negative 
time means before the present. 

the chosen volume element. The gas motion in the Z-direction is governed by the force F z and 
the restoring force of the stellar disk F = -A 2 Z. It can be expressed by 

sMMt-Q) r' h cos(Mt-Q) C f ~ u 

Z(t) = Zo(t e )+ — — — F z cos{At)cIt — — F z sin{Af)dT (11) 

^ Jo & Jo 

, (see Eq. 5.713 of Chandrasekhar 1943 1. As initial conditions of the disk in Z we put Za = Zq — 

at the time origin f e = -570 Myr. Ignoring the effects of the (X, Y) component of the force 

associated with the mass accretion, and assuming circular rotation, we have 6(t) = 6 (t e ) + ^-(f - 

fe). 

The flow of HVCs moves parallel to the X-axis. Therefore, we can think of the interaction of 
the HVCs with the disk as an accretion front perpendicular to the X-axis, propagating in the X 
direction with a constant velocity v a . Then the position of the accretion front as a function of time 
is X a (t) = v a f + Xo(f e ) (see Fig. l9-a> . The volume element enters the front when X — R sin(6(t)) = 
X(f) a , a condition that allows us to determine the entry time t \ . The time at which the volume 
element leaves the rear zone of the accretion front, t^, is determined by solving X = X(f) a + AX a 
for t, where AX a is the thickness of the accretion front. The times t \ and tj and the duration of the 
process of mass accretion, f a = t^— 1\, depend on 9 and R of the volume element (see Fig. 110-aL At 
large galactocentric radii the time intervals of mass accretion f a in regions of the third and fourth 
Galactic quadrants are longer than those in the first and second quadrants, i.e the outer regions in 
the third and fourth Galactic quadrants are the most affected by the interaction with the HVCs. 
For each volume element, the force F z acts between the corresponding times t\ and ?2, vanishing 
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Fig. 9-b. Angles (filled and open circles) and velocities (crosses) of entry of the test particles into 
the Galactic layer versus the respective positions projected on the X-axis. The open circles refer 
to angles of entry from above the Galactic plane. 

outside this time interval. Hence, Eq. (1 1) becomes Z(6,R,f) = ^{cosA(t2+t e -t)-cosA{t\+t e -t)). 
This is a general expression for the Z deformation of the Galactic layer due to the flow of HVCs. 
At t=0 this equation can reproduce roughly the present configuration of the Galactic warp, if 
the period of vertical oscillation Jr? ~ 450 Myr (see Fig. 110-bl . This value seems reasonable, 
because at large galactocentric radii the restoring force should be small. Note that the model 
can explain the overall pattern of positive vertical displacements in the first and second Galactic 
quadrants, as well as the less pronounced displacements toward negative Z in the third and fourth 
quadrants (cf. our Fig. llO-bl with Figs. 57b and 60b of Burton 1991 ). 

The first explanation proposed for the Galactic warp was that the Galaxy is passing through 
a continuous intergalactic medium (Kahn & Woltjer 1959 1. In contrast, we propose an interac- 
tion of the Galactic disk with a transient circumgalactic flow of HVCs produced by a single 
event of mass transfer from the Magellanic Clouds. Our model shows that the flow of HVCs 
is able to produce the formation of the Galactic warp. Future extensions of the model should 
contemplate consequences of the process of mass accretion such as large deviations from the 
circular rotation of the gaseous disk and the sudden increment of the rate of star formation and 
of the input of energy in the interstellar medium. The HVCs have injected energy at a rate of 
~ 10 50 ergs kpc~ 2 Myr -1 into the Galactic disk during the last 400 Myr. A question we should 
address, among others, should be the role played by this important source of energy in the pro- 
duction and maintenance of the diverse phases of the interstellar medium (Kulkarni & Heiles 
1988 Elmegreen 1991; Heiles & Troland 2003 1. Another interesting question is whether there 
exists any genetic or dynamic link between the Galactic warp and a ring of stars in the plane of 
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Fig. 10-a. The time t \ at which the head of the accretion front crosses a disk region, and the time 
?2 at which the back of the accretion front leaves the region, as a function of the present azimuth 
for the regions at R = 22 kpc (dashed line) and R = 24 kpc (solid line). The duration of the transit 
of the accretion front, t & = ti — 1\, is also represented as a function of the present azimuth. 

the Galaxy, at a Galactocentric radius of 18-20 kpc, which may completely encircle the Galaxy 
(Ibata et al. l2003l Yanny et al.[2003 ). 

6. Conclusions 

We have shown that the majority of the HVCs have a common origin linked to the dynamic 
evolution of the Magellanic Clouds. During a long time these two small satellite galaxies have 
likely behaved as a binary system under the effects of the tidal perturbations of the Galaxy. The 
two Clouds have probably approached each other several times in the past, until the Clouds 
collided catastrophically about 570 Myr ago. We assumed that the collision triggered a burst 
of stellar formation in the Clouds. The superwind driven by the starburst blew the HVCs as 
magnetized clouds into the massive halo of the Galaxy. As a consequence, the Milky Way is 
likely surrounded by a mist of highly ionized gas at high velocity, the HVCs being the denser 
concentrations. 

The catastrophic loss of matter was of a significant amount. Probably it permitted the SMC 
to reach escape velocity, disrupting the binary system. According to the results of our model the 
amounts of matter and mechanical energy liberated in the form of HVCs were of the order of 
3.1 x 10 9 Mq and 1.8 x 10 57 ergs, respectively. In the numerical simulation of the process, we 
assumed that the HVCs were ejected from the Clouds 570 Myr ago in all directions with peculiar 
velocities larger than the escape velocity of 160kms _I , and following a Gaussian distribution 
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Fig. 10-b. Vertical Z-distances of the midplane of the Galactic layer from the plane b = at the 
Galactocentric distances R = 22 kpc (dashed line) and R=2A kpc (solid line) as a function of the 
azimuth of the region, according to the simulation of the interaction between the accretion front 
due to the flow of HVCs and the Galactic gas layer. 

law. Because of these initial expansion velocities, the cloud of HVCs, or metacloud, expanded 
from the center of collision of the Magellanic Clouds, reaching its present size of ~ 300 kpc, and 
filling a great part of the halo volume of the Galaxy. As the initial mean velocity of the HVCs 
was that of the Clouds at the time of the ejections of HVCs, the centroid of the group of HVCs 
(metacloud) describes an orbit close to the orbit of the Magellanic Clouds (i.e. approximately in 
the direction I = 90° - 270°). To an observer at the position of the Sun, this gives the impression 
of a flow of HVCs coming from t = 90°. As a consequence of the initial velocity distribution of 
the HVCs, a larger number of HVCs tended to concentrate towards the Magellanic orbit behind 
the present position of the Clouds, whereas other HVCs moved ahead of them, forming the 
Magellanic Stream. 

The passage of the HVC flow through the Galactic disk had transcendent consequences for 
the evolution of the gaseous layer of the Galaxy. The outer gas layer of the Galaxy was consid- 
erably displaced in the Z-direction from its equilibrium position. The Galactic disk accumulated 
mass coming in the form of HVCs at an average rate of 6OOM0kpc~ 2 Myr~' over a period of 
200 Myr. Due to the time dependence of the mass accretion upon the azimuth of the region, the 
effects on the first and second Galactic quadrants at large galactocentric radii were different from 
those on the third and fourth ones. These are perhaps some of the clues to understand the gene- 
sis of the Galactic warp. The collisions of the HVCs with the Galactic disk might have exerted 
notable influences on the large-scale morphology and energy of the interstellar medium. 
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